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STATICS  AND  STABILITY  OF  THIN  WALLED  ELASTIC  BEAMS 


ABSTRACT 

Formulation  of  fundamental  equation*  of  elastic  equilibrium  of  thin 
vailed  beaus  subject  to  general  loads  and  dislocations  starting  only 
fro*  the  hypothesis  of  non  deformed  transverse  cross  sections. 
Formulation  of  the  fundamental  equations  of  dynsmic  stability  of  thin 
vailed  beams  subject  to  general  conservative  loads  and  dislocations 
by  use  of  a  systematic  geometrical  approach. 
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1)  STATICS 


1,1)  Introduction 

Th«  theory  of  elastic  equilibrium  of  a  cylinder  subject  to 
loads  applied  at  the  bases  and  represented  by  a  general  ays tea  of 
balanced  forces,  has  been  accurately  and  completely  developed  by  St. 

Venant  MOO  with  the  traditional  hypotheses  of  om^.enelty,  isotropy 
and  linear  elasticity. 

This  study  represents  the  background  of  the  so  called  "techaical 
theory  of  the  besas"  which  applies  with  approziaation  the  results  obtained 
by  St.  Venant  to  all  the  real  cases  concerning  the  elastic  equilibrium 
of  cylinders  subject  to  any  type  of  loads  and  constraints. 

Such  application  is  founded  in  a  classic  postulate  carrying 
St.  Veuant's  nans  and  is  synthetically  expressed  by  the  following 
principle:  "if  a  system  of  balanced  forces  acts  on  a  Halted  area  S' 
of  the  surface  S  of  a  body,  its  effects  danp  nut  as  they  leave  S*  end 
actually  UiSappear  at  distance  jp  depending  upon  the  shape  and  the  sis* 
of  s*': 

Such  postulate  permits  to  determine  stresses  and  displacements 
Having  knowledge  only  of  six  classic  stress  characteristics  connected 
with  the  constraints  and  Ji^ada  applied  to  the  body,  the  areas  close  to 
constraints  or  concentrated  loads  excluded. 

However,  soma  conditions  are  indispensable;  of  which  the  no at 
important  are: 

1)  the  cross  section  dimensions  must  be  coaparable; 

2)  the  body's  length  must  be  auch  greater  than  the  above  mentioned 
cross  dimensions. 

It  is  the  classic  case  of  solid  section  beams  for  which  the 
technical  theory  has  a  good  correspondence  with  reality. 

The  same  thing  does  not  ajiply  to  thin  walled  beams.  In  fact, 
auch  structures  are  characterized  by  three  diaenalonc,  anyone  of  which 
ia  negligible  if  compared  to  the  next  one: 

a)  thickness  of  the  wall 

b)  average  distension  of  the  cross  section 

c)  length 

For  this  type  of  structures,  which  are  always  more  widely  used 
hgr  technical  practice,  it  has  been  necessary  to  generalize  the  results 


obtained  by  St.  Venant,  specifically  as  far  as  torsional  stresses  are 
concerned;  a  new  theory  has  bsen  expressed  justifying,  with  approxi¬ 
mation,  the  discrepancies  between  technical  theory  and  te2t  controls. 

This  new  theory  known  as  "the  theory  of  sectorial  areas', 
developed  by  VlasovQs^OO  snc*  Timoshenko  [5]  DO  for  bee as  of  open 
croas  section,  has  been  laver  general itod  by  Karmen-Christensen  fVJ 
for  teams  of  general  cross  section . 

Vlasov  fjy  [4]  ,  Wagner  ,  Kappua  £9],  Goodier  |\GJ ,  etc. 
applied  this  theory  to  the  problem  oC.  elastic  equilibrium  stability 
and  their  results  have  been  confirmed  by  t*3t  controls. 

Nevertheless,  as  It  has  been  noticed  by  Karman-Vei-Zang- 
Chien  M.  vhe  sectorial  area  theory  is  only  the  fir.it  term  of  a 
repetition  procedure  the  validity  of  which  is  in  certain  case'*  doubtful. 

Such  theory,  !r«  fact,  basically  consists  ia  dividing  the 
shear  fto*  produced  by  the  twisting  moment  into  t¥©  peris:  tbs  priEary 
snear  flow  typical.  e£  3t.  Venant ’a  theory,  sod  the  saeoadery  shear 
fio*  associated  with  tae  aomal  ctressss  caused  by  the  son-i-nifors 
warp-Lag  of  eroaa  auctions  duo  to  the  primary  flow. 

In  fact,  the  sectorial  areas  theory  relapses  into  ths  classic 
solution  of  St,  Veasnt  whan  warping  is  ©easts at  in  the  length  ox  th@ 
bssss.  However,  this  theory  neglects  the  warping  caused  by  the  secondary 
sheer  flow  which  sosstisses  ess  be  ssrs  conspicuous  than  the  primary  one, 
sad^  coEieque&ily ,  t undasantally  changes  the  static  condition;  further¬ 
more,  said  theory,  oven  improving  considerably  the  correctness  of  cal¬ 
culation  of  stresses  and  deformations  inside  the  body,  cannot  be  applied 
ia  the  areas  which  srs  close  to  constraints  or  concentrated  loads. 

A*  e  conclusion,  ws  can  say  that  the  "sectorial  areas  theory" 
is  for  the  thin  walled  beams  what  th®  "technical  theory"  is  for  the 
solid  section  beams;  In  other  words,  the  limitations  of  both  can  be 
considered  identical. 

Therefore,  in  thin  Note  we  want  to  re-examine  from  the  origin 
the  p rob lest  of  a lagtlcequi librium  of  thin  wal lwd-be**B  ns  abject  to  very 
general  loads  and  dislocations,  making  use  of  a  very  general  method. 

In  fact,  the  correct  solution  of  the  problems  permits  to 
eliminate  the  limitations  related  to  St.  Venant’s  postulate  and  to 
determine  exactly  some  problems  of  considerable  Interest  for  the  theory 


as  veil  e.i  for  tho  practice,  as: 

1)  the  calculation  of  stresses  In  the  areas  close  to  concentrated 
loads  and  external  constraints; 

2)  the  calculation  of  stresses  associated  wilth  general  loads  acting 
on  the  surface  of  the  bean; 

>.ne  calculation  of  stresses  associated  1th  general  dislocation, 
of  general  interest  for  the  study  of  thersic  or  plastic  actions. 

In  the  first  part  of  this  study  the  prob.lem  of  elastic  equi¬ 
librium  of  thin  vailed  beams  vill  be  considered  froa  a  general  viewpoint 
sad  basic  equations  and  boundary  conditions  vill  be  furnished;  then, 
above  santioned  problems  vill  be  studied  and  solved. 


1,2)  The  basic  hypothesis 

The  basic  hypothesis  on  vhich  we  found  our  study  is  the 
hypothesis  of  s  transversal^  Sndeformable  cross  section.  Such  hypo¬ 
thesis  which  appears  also  in  the  theory  of  sectorial  areas  and  In 
Herman's  study,  is  generally  acceptable  for  the  thin  vailed  beams 
because  of  shear  dlaphrams  used  fcr  structures  of  this  type  with 
the  purpose  of  avoiding  the  buckling  of  the  vail. 

Such  dlaphrams  are  usually  realized  by  means  of  thin  plates 
welded  to  the  vail,  in  order  to  avoid  deformations  of  tho  cross  section. 

Nevertheless,  being  such  pistes  very  thin,  ve  can  Imagine 
them  having  no  resistance  to  varplng  outside  their  plane,  and,  conse¬ 
quently,  leaving  the  beam  cross  section  free  to  varp. 

Therefore,  In  this  study  ve  vill  consider  the  profile  as 
uniformly  stiffened  along  Its  vbole  length,  that  is,  ve  will  consider 
every  section  a«  keeping  unchanged  its  shape  during  the  displacement 
associated  with  general  loads  conditions. 


1,3)  Kinematic  relations 


VI th  reference  to  the 
thickness  t  and  a  general  cross 
0  the  shear  center  of  the  cross 


profile  shovn  in  fig.  1,  haring  a  constant 
section,  ve  denote  G  the  centroid  and 
section. 


3 


■» 


We  refer  the  points  of  the  surface  to  the  orthogonal  right-hand 
tern  Gxyz ,  of  which  axes  x  and  y  coincide  with  the  principal  inertia 
axes  of  the  cross  section  and  axis  z  is  perpendicular  and  passes  through 
the  centroid  G.  Furthermore,  we  refer  the  bean  surface  to  the  two  groups 
of  orthogonal  lines  foraed  by  directrices  and  generatrices  of  the  cylin¬ 
drical  surface,  choosing  n  normal  to  the  surface  in  a  point  P(s,z),  and 
e  and  z  such  that  the  directions  tern  (n,  s,  z)  is  right  and  can  be  super¬ 
imposed  on  fixed  tern  Gxyz  with  a  rigid  motion. 

Being  ^  the  displacement  of  point  P,  we  denote: 

u  «  u  (x  y  z), 

v  -  v  (x  y  z)  ,  (1 

w  -  w  (x  y  z) , 

the  components  of  such  displacement  on  the  axes  of  fixed  tern  xyz,  end 
we  demote: 


The  basic  hypothesis  permits  to  determine  the  displacement 
in  the  plane  z  y  of  every  point  of  the  cross  section  with  only  three 
parameters  only  depending  upon  abscissa  z.  In  fact,  denoting: 

“o  =  “o  <z> 

vo  x  vo  (1’7> 

f*  -  fm  <*> 

the  displacement  component  on  x  and  y  of  the  shear  center  0  and  the 
section  rotation  around  the  stuar  center,  the  first  two  equations  (1,1) 
can  be  written  as  follows: 
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Cl,  8) 


u(*.f.  *)-  *•(*)- 
f-fc f.zj-  KfcJj-  </„(£■  x.J 

where  (xc,  yD)  are  the  coordinates  of  the  shear  center  0  (fig.  2). 
Therefore,  using  equations  (1,6),  we  have: 


Fig.  2 

We  observe  that  the  quantities: 

r  *  "  xo)  S  '  (y  '  yo}  S  ' 

”  (i,io) 

p  - <x  •  V  s +  (y  -  *>>  S  * 

are  the  components  on  the  a^es  a  and  n  of  vector  R  =  OP,  therefore 
equations  (1  2)  can  be  finally  expressed  as  follows: 

gl£ +  cf.t  (i,ii) 

0*3  ot  3 
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Consequently,  the  ration  of  every  point  of  the  beam  is  expressed 
by  the  following  four  functions: 

Uo  *  u0  • 

vo  *  vo  (z)  * 

w  =  w  (n.s.z)  , 

and  the  latter  can  be  considered,  with  a  good  approximation,  independent 
of  n,  in  consideration  of  the  smallness  of  thickness  t,  and  can  be 
written: 

w  *  w  (s,z)  .  (1,13) 


1 ,4)  Elasticity  relations 


If  we  neglect  the  normal  stress  CJj  ,  we  can  express  as  follows 
the  relations  between  the  stresses  components  and  the  unit  strains  In 
the  thin  wall  surface: 


0 


(1,14) 


^ — '  /Va 


In  a  rare  general  case  the  strain  components  will  be  expressed 
by  the  following  relations: 


mr  4^ 


(1,15) 


where  &  Is  the  elastic  strain  and  £  the  strain  due  to  a  general 

dislocation  system,  as  a  thermic,  plastic  system  etc. 

So  equations  (1,14)  can  be  written: 

/sr_  -J.  if  I  JF  _  /, r-  f  ,  U 


^  a  *  f*rJ 


/ms - 

'  2{U*J  * 


(1,16) 
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Furthermore ,  the  basic  hypothesis  permits  to  reduce  the 


unknowns;  in  fact,  since  we  must  have: 

=  0 


(1,17) 


for  the  cross  indeformability  of  the  section,  the  normal  stress  can 


be  expressed: 


while  the  shearing  strain  /-•  can  be  expressed  as  follows: 

* 

where  ,  ijy  ,  and  are  general  functions  of  z. 

The  basic  unknowns,  expressed  as  special  stress  components, 
can  be  reduced  to  the  following  two  functions: 

~  *•) 


(1,18) 


(1,19) 


where : 


a  »  m 


Taking  into  account  the  classic  relations: 


(1,20) 


(1,21) 


(1,22) 


equations  (1,20)  ctn  be  written  as  follows  for  (1,6): 


and  they  express  the  general  elasticity  relations  of  thin  walled  beams. 

Equations  (1,23)  represent  the  values^ of  normal  stresses  £5^ 
and  shear  stresses  corresponding  to  the  middle  fiber  of  the  wall 

forming  the  profile. 

In  reality  such  stresses  vary  along  thickness  t  of  tho  wall, 
but  actually  they  can  be  considered  constant  because  of  the  thickness 
smallness.  However,  if  the  profile  has  open  cross  section,  it  is  neces¬ 
sary  to  consider,  together  with  the  stresses  (1,23),  the  shearing  stres¬ 
ses  linearly  variable  along  the  thickness  and  vanishing  in  correspondence 
with  tho  middle  fiber  associated  with  the  twist  of  the  wall  caused  by 
external  torque. 


3ucb  stresses,  classic  of  St.  Venant’s  study,  can  be  expressed, 


with  good  approximation,  as  follows: 


zfi.J 


(1,24) 


being  n  ihs  distance  becween  the  fiber  and  the  Biddle  surface;  in  fact, 
said  stresses  are  the  only  ones  which  develop  for  a  constant  twist  of 
the  beam  and,  consequently,  allow  the  beaa  to  balance  the  external 


torque. 


In  fact,  as  a  result  of  (1,24)  we  obtain  a  twisting  moment 


M  ,  having  the  well  known  expression: 

Z 


(1,25) 


being  JT*  the  torsional  rigidity  which,  in  case  of  open  sections  of 


constant  thickness  t  ,  is  written: 


X-  ml3 

3 


(1.26) 


where  m  is  the  length  of  the  middle  line,  and  in  the  case  of  cross  section 
consisting  of  several  portions  of  different  thickness  tA,  is: 

f*  5 

JE.  rrt;  A-  (1,27) 


If  the  profile  ha^  a  close  section  (box  or  multicell  beam), 
stresses  (1,24)  are  no  more  necessary  to  give  torsional  rigidity  to  the 
beam.  In  fact,  also  in  case  of  constant  twist,  the  external  moment  is 
almost  completely  absorbed  by  a  flow  of  shear  stresses  constant  along 
the  thickness;  and,  compared  with  such  stresses,  the  contribution  given 
by  equations  (1,24)  is  quite  unimportant. 

Therefore,  in  these  cases,  stresses  (1,23)  are  sufficient  to 
balance  any  external  action  and,  consequently,  are  the  only  stresses 
which  are  considered  acting  on  the  wall. 


1,5)  Ilqullibrlua  equations 

With  reference  to  the  wall  element  ds  dz  inside  point  P(s,z) 
of  the  middle  surface,  the  equilibrium  equations  to  be  imposed  coincide 
with  the  three  equilibrium  conditions  relative  to  the  displacement  along 
axes  n,  s,  z.  The  first  two,  concerning  the  equilibrium  along  normal  n 
and  tangent  s,  become  unessential  because  of  the  hypothesis  on  the  inde- 
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fornability  of  cross  section  of  the  bean.  In  fact,  in  such  directions 
the  equilibrium  is  guaranteed  by  the  mutual  actions  of  the  stiffeners 
on  the  wall  which  can  be  so  calculated. 

Therefores  if  we  denote  p2,  px,  py,  (fig.  3) .respectively , 
the  load  acting  in  the  direction  of  axis  z  on  the  wall  element  ds  dz; 
the  loads  acting  in  the  direction  of  axes  x  and  y  and  the  twisting 
mcoeat  on  an  element  of  the  beam  having  length  dz;  the  equilibrium 
equations  are  written: 

V-  ,  JE,  «  o 

+  o 

pi-O  (1,28) 

a/x  * 


•fa* 


c/r: 


+  Mm- 


o 


Fig.  3 


being  T^,  T  ,  Mz  the  resultants  of  internal  stresses  in  the 

direction  of  axes  x  and  y  and  the  resultant  moment  in  regards  of  shear 
centers  0. 
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These  letters  can,  therefore,  be  expressed  as  follows: 


(1,29) 


being  the  Internal  moment  expressed  by  equation  (1,25)  and  associated 
to  stresses  which  will  be  taken  into  account  only  in  the  case  of 

open  sections. 

The  last  equation  (1,29)  can  have  the  sane  form  for  open 
sections  as  well  as  for  box  or  multicell  sections,  by  introducing  a 
warping  function  associated  with  constant  twist. 

Such  function,  which  we  denote  COm  ,  represents  the  axial 
warping  function  w(s)  of  the  points  of  the  wall  middle  line  when  subject 
to  a  constant  torque  having  unitary  negative  gradient 

in  the  case  of  open  section  beams,  such  function  is  obtained 
by  observing  that,  since,  in  accordance  with  St.  Venant’s  solution, 

2? MB  equals  0  in  correspondence  with  the  middle  line,  the  second 

equation  (1,23),  having: 

Mffa  Sjw  (*J 

O  (1,3°) 

afr  * 


gives : 

0 

(1,31) 

On  the  contrary,  in  the  case  of  close  or  multicell  sections 
(fig.  4),  such  function  can  be  obtained  by  considering  that,  since  T^xe 
coincides  with  the  flow  of  stresses  resulting  from  known  solution  of 
Bredt-St,  Venant: 

2ws»  &&£ 


(1,32) 


being  f  the  flow  constant  (*) ,  the  second  equation  (1,23),  in  conside¬ 
ration  of  (1,30)  and  (1,32),  gives: 

o 

t 


(1,33) 


Fig.  4 


(*)  We  must  remember  that  flov  constants  t,  for  every  element  of  multicell 
section,  can  be  obtained  with  the  partial  flow  and  relative  to 
meshes  1  and  k  having  such  elements  in  conaon.  The  partial  flov  constants 
1 1  can  finally  be  obtained  from  mo  nod  romi c  condition  of  and,  con¬ 

sequently,  of  6t?m  which  Imposes  for  every  circuit  the  following  relation: 


jfta/s  — j.£  $/&  .«•  O 


from  which,  denoting  8,  the  area  enclosed  by  circuit  i,  we  obtain: 

X 
rmt  ' 


y&t'-  y-  2E'f0cf;ie  »  O 


(1 ,33) 


(1,34) 


where  Off  represents  the  geometric  circuitation: 

at'/mm  o/x 

relative  to  the  whole  circuit  i  and  O^/tc  represents  the  partial  geome¬ 
trical  circuitation  of  the  element  in  common  to  meshes  i  and  k.  Eqs.  (1,34) 
represent  a  system  which  is  linear  for  unknowns  f^  and  of  simple  solu¬ 
tion.  In  view  of  the  above  it  is  easy  to  obtain  constants  f. 
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Therefore,  from  equations  (1,31)  and  (1,33),  with  a  simple 
quadrature  procedure,  we  can  obtain,  neglecting  an  arbitrary  constant, 
the  expression  of  function  OtJm  .  In  general  the  constant  is  elimi¬ 


nated  with  the  auxiliary  condition: 


o 


Thus,  equations  (1,23)  can  be  written  as  follows: 


(1,35) 


o/cc/m 


/J 

<3/ Si  . 


I o/jr  / 


(1.36) 


for  open  sections,  and: 


for  close  or  aulticell  sections. 

Denoting  ftl  the  nuaber  of  close  aeshes  of  cross  section,  from 
equation*.  _X1^23)  we  have: 


(1,38) 


(1,39) 


(1 ,40) 


and,  Introducing  the  notation: 


(1.41) 


•^uatlM  (1,37)  can  be  written  as  follows: 

K**  fir9m  oAv*o//4  + 

/  o/s  l o/z  * 

4 

and  it  appears  identical  to  tbe  equation  already  obtained  for  open 
sections  and  expressed  by  (1,36), 

Therefore,  without  considering  the  type  of  bean  cross  section, 
equations  (1,29)  can  be  written  as  follows: 


/ 


(1,42) 


resulting  connected  to  the  cross  section  geometry  by  the  three  basic 
functions : 

x  *  x  (s)  y  =  y  (■)  cOg  =C&J.s) 

Por  these  functions,  we  nust  receober  that,  since  we  chose 
axes  x  and  y  as  aain  inertia  axes  and  the  center  or  rotation  0  as  shear 
center  of  the  section,  we  will  always  have  the  basic  relations: 


jxya/J=jxa>c 


& 


(y  CO*  o 


(1,43) 


In  fact,  we  can  ootaln  the  coordinates  xQ  and  yQ  of  ceni  u  -t 
rotation  0  by  laposing  the  last  two  equations  (1,43),  or  by  using 
Jouravsky's  procedure  for  close  or  cellular  sections  £12}  . 
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1,6)  Basic  equations  of  elastic  equilibria*  of  thin  walled  beam  having 
continuous  directrix  and  constant  thickness 


Ve  can  no*  obtain  the  basic  equations  of  elastic  equilibrium  of 

thin  vailed  beams,  by  changing  the  indefinite  equilibrium  equations  (1,28) 

into  terms  of  displacement.  For  the  moment,  since  *e  consider  the  body 

free  in  the  space  snd  subject  to  a  system  of  balanced  forces,  we  know 

the  three  transversal  characteristics  T_  (z) ,  T_  (z)  ,  M_  (z)  ,  and  »e 

*  y  z 


can  simplify  equations  (1,28)  as  follows: 

£5+  —ii  +  e±  m.O 

<*&  tr 

•//«.  7r 

T* 

+  <PX^z(„  &mJm  //e 


(1,44) 


The  first  of  these  equations  expresses  the  equilibrium  in  the 
direction  z,  and  the  three  other  ones  express  the  identity  between  the 
resultants  of  internal  shear  stresses  and  the  corresponding  stress 
characteristics;  such  equations  can  be  expressed  for  the  displacement 
parameters  (1,12),  taking  into  account  the  elasticity  relations  (1,23); 
in  fact,  we  have: 


f  * 


alter  Introducing  the  notations: 

~  */•*  ~ f-jr~  j?*9 


(♦♦)  (1,46) 


Therefore,  equations  (1,45)  are  the  requested  elastic  equilibrium 
equations  of  thin  walled  beams  subject  to  loads  and  dislocations.  Such 
system  can  be  simplified  by  drawing  from  last  three  equations  the  functions 

in  function  of  the  stress  characteristics  Tx,  Ty, 
and  axial  displacement  HZ' For  t  ils  purpose,  denoting  D  the  determi¬ 
nant  of  symmetrical  matrix  coefficient- 

(U,  •/(!  • 

J"  Jt‘  ‘‘"I  (1,47) 

«/|»  </»•/ 


(1,47) 


and  denoting  the  complementary  matrix  of  element  d^^,  we  obtain  from 
las:  three  equations  of  system  (1,45)  the  following  relations: 


(**)  The  last  two  equations  (1,4)  are  directly  verified  for  open  sections 
with  equation  (l,3f:)  and  for  close  or  cellular  sections  with  equation  (1,33) 
and  with  the  following  relations: 

° 

which  ensue  from  the  equilibrium  condition  in  the  x  and  y  direction. 


+  -ff-  Zb*  M* J 


If  we  operate  in  the  second  part  of  equations  (1,48)  the  following 
linear  changes: 


v„yt*u 

y fPt  x/i&J /-  pgty /  sj  j.  Pta  co'a  C&J j 
£*  [ sj- £ ^ Pa  */&;+  P„ yfrj-t  Ppf  CCJs  <*>/ 


(***)<!, 49> 


(***)  Equations  (1,49)  require  the  following  invorsed  relations  to  be  true; 

x(sj*  o/„y^ 


*•/»  1/„ y'xj-s- 

J  (1,49)’ 

COtC&jrro/a  Y(sj-J-  y/fej-/-  s/at 

furthermore,  it  is  easy  to  verify  that  the  six  functions: 


have  the  following  properties: 


*£/*//  f/J?: 

o/s  o/s  *  ’&7/S 


) 


fjy  </$• 

Jd&  3s 

J* 


slim r  0 


3S  (1,49)" 


which  can  be  controlled  taking  into  account  the  determinants  properties 
and  the  relation: 
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(1.501 


+  T>,'  //./ 

/ /3, T,  +  P..7j+  ?.. M.j 

fp*  ^  ^2.  3^"’* 

we  can  simplify  as  follows: 

«r  £<ts  a/s  e 


a,;rt 


5^-r-  'y  ffljgj.  ifl 

ofr  J/*  <*''•*  & 

Therefore,  the  basic  equation  of  thin  walled  beans  is  obtained 
by  substituting  (1,51)  in  first  equation  (1,45)  and  observing  that  for 
equations  (1,31)  or  (1,33)  we  always  have: 


o/cL>c  a.  &/g 

a/ *  a/s  (i.s2) 


in  consideration  of  the  above  and  taking  into  account  equations 
(1,49)  and  tl,50>,  we  obtain  the  following  integral  differential  linear 


which  is  of  basic  importance  for  the  study  o'  thin  walled  beams  of 


tr&nsversally  indeformable  section  subject  to  general  forces  and  dislo¬ 
cations.  In  view  of  future  applications,  it  is  therefore  adv  able  to 
express  the  elasticity  relations  (1,23)  by  the  displacement  ax. pI  compo¬ 
nent  This  can  be  done  simply  by  taking  into  account  eqi  itions 
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(1,51);  therefore  we  have: 


(1,54) 


which  are  the  final  expressions  cf  elasticity  relations  for  thin  walled 
beau.  Equation  (1,53)  Bust  furnish  solutions  satisfying  the  boundary 
conditions  on  the  bases  (z  =  0  and  z  *  1)  and  the  transversal  conditions 
depending  upon  the  ah  ye  of  the  bean  section  described  in  following 
paragraph. 


1,7)  Boundary  conditions  connected  with  basic  equation 


Ve  divide  the  boundary  conditions  into  longitudinal  conditions, 
regarding  the  external  bases  z  *  0  and  z  =  1,  and  transversal  conditions. 

In  case  of  longitudinal  conditions  we  notice  that,  if  ve  con¬ 
sider  a  body  fiee  and  subject  to  a  systes  of  balanced  forces,  said  con¬ 
ditions  rill  necessarily  iapose  the  equality,  in  every  point,  between 
external  actions  and  acting  respectively  on  bases 

z  =  0  and  z  =  1,  and  corresponding  normal  stresses  and  &Z s) ; 

therefore,  they  are  as  follows: 


crM[lts)=*  jz>49 


(1,55) 


Equations  (1,55,  expressed  with  equations  (J 4)  for  displacement  give: 


which  represent  the  two  necossary  longitudinal  conditions  to  be  associated 

wi *■  basic  equation  (1,53).  We  notice  that  on  bases  z  «  0  and  z  =  1 

the  iaentity  in  every  point  between  external  actions  and  internal  stresses 
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concerns  only  normal  stresses  and  not  sheer  stresses  for  which  equations 
(1,44)  guarantee  global  identity  referred  to  resulting  actions  (forces 
and  oosent) . 

As  far  as  the  end  bases  are  concerned,  the  difference  in 
every  point  between  external  actions  and  internal  stresses 

2^,*  is  entirely  absorbed  by  two  existing  stiffeners  and,  consequently, 
does  not  cause  any  additional  deformations  or  stresses  not  even  In  the 
areas  very  close  to  the  two  bases. 

Equations  (1,56)  are  therefore  the  only  longitudinal  conditions 
concerning  the  extreme  bases. 

A  different  procedure  is  required  for  transversal  conditions, 
since  they  depend  upon  the  type  of  the  cross  section.  Therefore,  we 
will  consider  '  hem  case  by  case  in  regards  to  the  sh  _■  of  the  cross 
section  direct* 

a)  Open  sections  having  continuous  directrix 

Ve  consider  as  continuous  directrix  a  curve  having  functions 
x(s),  y(s)  ana  co  (s)  continuous  vp  to  the  second  derivatives;  such 
sections  (fig.  5)  cannot  have  aor >  than  two  generatrices  and  we  denote 
a j  and  s2  respectively  their  curvilinear  abscissa. 


Fig-i — 5  —  - 

If  we  denote  <?”.)  and  P2(z)  the  tangential  loads  eventually 
acting  on  such  generatrices,  and  t  the  constant  thickness  of  the  wall, 
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we  can  write  the  transversal  conditions  as  follows: 


7%*J 

S' 

2i  .(S.,z)= 

*- 

o 


( 1 , 57)  a 


Taking  into  account  equations  (1,54)  and  denoting  L(w)  the 


term: 

c?s 


-  v  f */J 

Jc?s.  a/s 

'A 

equations  (1,57)  can  finally  be  written  as  follows: 


(1,58) 


‘■(“U-  ^r‘V 


(1 ,59)a 


b)  Close  sections  having  continuous  directrix 

In  addition  to  what  stated  in  paragraph  a)  above,  concerning 
the  definition  of  continuous  directrix,  for  these  sections  (fig.  6) 
the  transversal  conditions  will  be  expressed  as  continuity  conditions 
for  functions  sj  and  Sj  (being  t  constant)  in  the 

limited  field  of  curvilinear  abscissa  s.  Such  conditions,  reflecting 
the  double  aspect  of  geometrical  compatibility  and  equilibrium,  will 
be  expressed  as  follows: 


which,  because  of  equations  (1,54)  and  the  hypothesis  of  continuous 
coordinate  functions,  become: 


(1 , 57) b 


(1 , 59)b 
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Fig.  6 


From  equations  (1,59)  in  the  form  a  or  b,  in  accordance  with 
the  type  of  cross  section  associated  with  the  longitudinal  conditions 
(1,56),  we  can  obtain  uni vocally  the  solution  of  basic  equation  (1,53). 

1 ,8)  Basic  equation  extended  to  thin  walled  beams  having  discontinuous 
directrix  and  discontinuous  constant  thickness 

In  reality  the  thin  walled  beams  are  nearly  always  formed  by 


more  than  one  element  (fig.  7) ,  everyone  of  which  can  be  considered 
as  an  elementary  beam  having  continuous  directrix  and  constant  thickness. 


Denoting  n  the  nunber  of  elements  forming  the  beam  and  i 
a  general  element,  the  displacement  parameters  indicating  the  motion  of 
every  point  of  the  beam  cross  section  will  be  the  n  +  3  functions: 

X.Cey 

r.(?) 

<1‘ 

<3, - f-J 

being  (si(  z)  the  axial  displacement  of  point  (sA ,  z)  of  the 
element  middle  surface. 

Then,  denoting  (s^ ,  (s  i>**fi  (s^)  ,  r^  (st)  the  functions 
typical  of  element  i,  and  (s^),  y^  (s^) ,  JZzi  (.t)  the  varied  expres¬ 
sions  : 

y,{sj-  J. '  fa,*. -(s;J+  X}s  cOt.(s,jJ 

where  Dik  are  always  the  complementary  matrices  of  elements  dlk  of 
determinant  (1,47)  which,  this  time,  we  express  as  follows: 


(1,12)  ' 


*~*@e>M* r 

<*/*, i  «■  j */Xj  &/J; 
a/,e  --T  r,  <s//z 

u. 

e/im  *  *,  o/d,- 


(1.49)' 


(1,46)' 
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the  equations  determining  displacements  LZ. 


sions  (1,50)  of  forces 


£ 


,  ,  </*>  , 

,  become : 


with  expres- 


__ 

o/z 


EEL 

‘V^/Vs/ 

V, 


,  J*XsO£-,/A-+ 

l***'*' 


(1,51)’ 


Therefore,  n  equations  determining  displacements  w;wil,l  be 
written  as  follows: 


where  t^  represents  the  constant  thickness  of  element  i;  pzi  (sA)  represents 
the  axial  load  acting  on  said  element  for  unit  of  surface;  and  «ST +  / 
represents  the  anelastic  strain  component  acting  on  the  same  element. 

Therefore  n  +  3  equations  formed  by  (1,51)'  and  (1,53)'  generally 
solve  the  problem  of  elastic  equilibrium  of  thin  walled  beams,  provided 
that  its  section  has  constant  discontinuous  thickness. 

In  order  to  solve  said  equations  we  must  find  the  longitudinal 
aad  transversal  boundary  conditions. 
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The  first  ones  express,  as  usual,  the  equilibrium  condition: 


(1,55)  ’ 


and,  reduced  in  terms  of  displacement,  give: 

fz&l '=»  —  ^ 

^ c) £  'mo  ^  < 

i (§fJ.T/  +  U 

while  the  transversal  conditions  will  concern  geometrical  compatibility 
and  equilibrium  conditions  corresponding  to  every  junction  point  of 
several  consecutive  elements. 

Denoting  k  the  number  of  elements  present  in  the  junction 
(fig.  8),  these  latters  will  be  written  as  follows: 


(1.56)' 


—  >&/  * j  J 

?;  r.*  zrjJ =-  7 


(1,57)' 


where  is  the  curvilinear  abscissa  of  the  Junction  in  relation  with 
element  i,  and  P1  (z)  is' the  eventual  external  tangential  action  acting 
on  the  junction  point  itself. 


fl  t, 


Fig.  8 


In  the  summation,  the  positive  signs  concern  the  elements 
having  curvilinear  abscissa  converging  in  the  junction  point,  and  the 
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negative  ones  concern  the  remaining  elements.  Therefore,  equations 
01 ,5-7) ' ,  written  in  terms  ef  displacement,  give: 


and  is  obvious  the  change,  if  the  end  of  element  i  is  free  rather  than 
connected  to  other  eleaents. 

1 , 9)  Conclusions 

Proa  the  study  perforated  it  anpears  clear  that  the  problem  of 
elastic  equilibrium  of  thin  vailed  beams,  considered  as  cylinders  having 
transversally  inde forma ble  profile,  is  acre  complicated  than  what  could 
be  expected  following  the  beams  technical  theory  or  the  store  recent  Vlasov's 
theory  of  sectorial  area. 

In  fact,  the  problem  can  strictly  be  expressed  by  an  integro-dif- 
ferentlal  equation  linear  to  the  partial  derivatives  in  unknown  function 
w  (z,  a),  which  physically  coincides  with  the  axial  component  of  points 
displacement  of  middle  fiber  of  the  wall. 

Such  equation  la  not  of  difficult  solution;  a  general  solution 
will  be  furnished  in  the  following  part  of  this  report,  showing  how  our 
solutions  are  similar  to  thoae  obtained  by  abo/e  mentioned  approximate 
theories  and  pointing  out  the  unavoidable  approximation  of  same  theories. 
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2)  STABILITY 


2,1)  Introduction 

T^MKjggigKntnl  static  behaviour  o f  a  thin  vailed  beam  In  regards 
of  St.  Ven^^^Pj^p^der  appears  conspicuous  when  we  study  Its  equilibrium 
stability the  classic  behaviour  of  the  beam  subject  to  combined 
bending  gB^Pompressl  ve  stress,  which  bends  In  a  main  inertia  plane  of 
the  sectlojSg^t  caso  of  unstability  due  to  axial  stress,  the  thin  vailed 
beam  of  bends  and  tvlsts  at  the  same  time  under  loads  such 

analler  ti^^d^T^correspODding  to  Eulero’s  formula. 

^^jgggggj^blllty  of  having  a  f lexlo-torsional  buckling  under 
axial  st£g£^riM*v**9£overed  when  thin  vailed  members  of  open  section 

rgt  time  In  designing  aeronautic  structures:  many 
,  Ostenfeld  Q. 3]  ,  F.  and  H.  Bleich  [h]  ,  and 
ted  the  lavs  governing  the  phenomenon.  Only  vith 
to  the  more  recent  studies  of  Timoshenko and 
her  with  the  works  of  Vlasov  [15]  and  Goldenveiser 
principles  have  been  established:  In  order  to 


axial  st£gf^riM*v**9£overed  when  thin  vailed  members  of  open  section 
were  use-""-- 
Authors 
Kappus  £ 

F.  and 
Goodier 

OJ  . 

determine  the  presence  of  bending  In  the  beam,  the  center  of  gravity 


had  to 


tuted  vith  the  center  of  torsion;  only  when  the  axis 


of  the  center  of  torsion  was  rectilinear  no  flexural  energy  was  present 
in  the  thin  vailed  members;  and,  furthermore,  the  warping  rigidity  Cj 
was  exactly  formulated. 

Of  great  Importance  are  the  studies  performed  by  Vlasov  [is] 
for  the  formulation  of  a  theory  concerning  the  unstabllity  of  the  thin 
vailed  beam  of  open  section  subject  to  normal,  bending  and  shearing 
stresses,  and  the  studies  performed  by  Krall  1.17],  who  obtains  the  stabi¬ 
lity  equations  by  using  the  variational  approach  with  the  Introduction 
of  the  twisting  moment  and  considering  vcrlous  cases  of  combined  unsta- 


blllty . 


The  constant  progress  of  technics  led  to  an  always  wider 


application  of  the  thin  vailed  beams;  this  structural  element  is  now 
present  In  most  civil  and  industrial,  naval,  aeronautic  and  space 


constructions 


Therefore,  the  study  of  equilibrium  stabi  ty  of  a  thin  walled 
beam  of  open  section  is  always  of  great  interest  and  new  problems  arise: 
as,  for  instance,  the  basic  one  concerning  the  influence  of  the  dislocation 


on  the  stability,  its  effect  and  the  effect  of  external  conservative  and 
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nonconsorvative  forces  on  the  dynamics,  etc. 

Thus,  we  want  to  examine  again  the  whole  system  of  elastic 
equilibrium  stability  of  the  thin  walled  beam  subject  to  general 
loading  and  dislocation  and  we  try  to  set  up  e  new  general  theory. 

The  study  of  such  beams,  as  conducted  in  the  first  part  of 
this  repor*  and  connected  with  researches  underlay,  the  results  of  which 
will  be  furnished  in  a  later  report,  confirm  the  validity  of  sectorial 
areas  theory  without  consideration  of  local  effects  connected  to  the 
presence  of  concentrated  forces,  holes,  etc. 

Thus,  Baking  use  of  Vlasov's  static  theory,  simple  and  suffi¬ 
ciently  correct  for  an  lnvest'gation  on  such  phenomena,  we  obtain,  in 
accordance  with  dynamic  method,  and  using  a  geometric  systeuatic  pro¬ 
cedure,  the  basic  differential  equations  governing  the  stability  problem 
of  the  beam  action  under  generally  distributed  conservati  e  forces  and 
dislocations.  Suci  equations  are  expressed  by  the  load:>  lirectly  applied 
and  the  stress  conpocents  corresponding  to  the  i»s.  c  cor  figuration  and 
includes  four  functions  characterizing  the  flexural,  toreional  inc 
extensional  oscillations  respectively.  The  extensional  oscillation  is 
often  neglected,  but  is  interesting  because  of  its  stabilizing  effects. 

The  systoB  of  forces  F°  (x  y  z)  acting  on  the  tain  walled  beam 

is  formed  by  distributed  forces  Q°  (z,  s) ;  Q°  (z,  s) ;  Q°„  (z,  s) ,  which 

x  y  z 

have  the  same  direction  of  axes  x,  y,  z,  and  are  functions  of  curvilinear 
abscissa  s  formed  by  the  center  line  of  the  cross  section.  Such  forces 
are  conservative  and  keep  their  direction  during  the  displacement  of 
the  points  at  which  they  are  applied  and  generally  originate  a  distri¬ 
bution  of  transversal  forces  p°  (z)  and  p°  (z) ,  axial  forces  p°  (z) , 

x  y  z 

bending  couples  m°  (z)  and  b°  (z) ,  twisting  couples  m°  (z) ,  and 
A  J  z 

bimoments  8°  (z) . 

The  dislocations  system  A,* (x  y  z)  causes  a  stress  condition 
which  can  be  annulled,  generally,  only  by  dividing  the  body  into  its 
elementary  particles  or,  more  simply,  by  cutting  it  into  a  finite  number 
of  planes.  The  Introduction  of  the  dislocations  system  (x  y  z)  will 
be  useful  later  for  the  study  of  the  unstabilizing  effects  caused  by 
residual  stresses,  non  uniform  thermic  field  or  prestressing  systems. 
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2 , 2)  General  remarks  on  approach  method 


Fig.  9  shows  the  axes  system  where  C  is  the  centroid;  x  and  y 

are  the  main  inertia  axes  of  cross  section;  z  is  the  centroid  axis. 

The  coordinates  of  shear  center  0  in  the  section  plane  are  x  and  y„ . 

o  o 


Fig.  9 


The  external  forces  are  generally  represented  by  components 

(z,  s) ,  Q°„  (z,  s) .  and  Q°_  (z,  s)  having  the  same  direction  of 
x  y  z 

axes  x,  y,  z  of  fixed  coordinates  system  Cxyz;  and  ara  general  functions 
of  curvilinear  abscissa  s  formed  by  the  center  line  of  thin  cross  section 
and  by  abscissa  z.  Such  forces  will  be  considered  as  conservative  forces 
and,  specifically,  as  keeping  unchanged  their  directions  determined  by 
fixed  axes  x,  y,  z  respectively.  The  loads  at  the  end  sections  are 
formed  by  a  distribution  of  general  forces  but  still  conservative  cor¬ 
responding  to  normal,  sheering,  bending,  twisting  and  warping  actions. 

The  coaction  state  due  to  dislocations  is  represented  by  normal  and 
shearing  stresses  in  every  cross  section  sel.’-balanced  if  the  external 
constraints  do  not  react. 

Therefore,  with  reference  to  a  general  cross  section  of  the 
body,  the  stress  state  will  be  represented  by  seven  stress  characteristics; 
bending  moments  (z)  and  (z) ;  twisting  moment  Mz  (z) ;  shearing  stresses 

Tx  (z)  and  Ty  (z) ;  bimoment  B  (z) ,  as  shown  in  fig.  10. 
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The  thin  beam  motion  will  be  formed  by: 

a)  a  system  of  displacements,  typical  jf  a  flexural  oscillation, 
by  which  the  axis  of  shear  centers  0  bends  in  the  planes  xz  and  yz, 
and  the  cross  sections  of  the  be  a  transfer  along  their  planes  and 
rotate  around  axes  x  and  y; 

b)  a  system  of  displacements,  typical  of  a  toisional  oscillation, 

by  which  the  cross  sections  transfer  along  their  planes,  rotating  around 
the  shear  centers  axis  (which  remains  rectilinear)  and  warp  because  of 
the  sectorial  areas; 

c)  a  system  of  displacements,  typical  of  an  extensional  oscillation, 
by  which  the  cross  sections  transfer  in  parallel  with  themselves  along 
the  direction  z  of  fixed  system  Cxyz . 

The  new  actions  developing  along  the  direction  z  on  the  element 
dA  dz  will  be  calculated  by  determining: 

1)  the  transversal  and  axial  elementary  forces  due  to  the  change  of 
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direction  of  stresses  f  4 A  ,  r  qA  .  9  d 'A  ,  folio®* ng  the 

**•  */  * 

fibers  buckling; 

2)  the  elementary  couples,  which  we  call  "turnover"  couples,  causing 
the  rotation  of  element  dA  dz  around  fixed  axes  x  and  y  and  z  rnd  due 

to  the  components  along  axes  x,  y,  z  of  fixed  systass  of  elementary  forces 

t-  6A  ,  v  dA  i  (T  dA  >  acting  on  the  buckled  body; 

ay  z 

3)  the  elementary  couples,  which  we  call  "displacement”  couples,  due 

to  the  fact  that  forces  'f  dA  ,  T*  dA  ,  dA  .  acting  on  the 

^  ^ 

two  sides  dA  of  element  dA  dz  and  the  surface  forces  Q  (zi)d- j,  Q  (x-jJ  d-i  , 
Q  z  f  during  the  buckling,  assume  a  different  position  in  regards 

of  fixed  referenced  system. 

Further,  we  calculate  the  actions  which,  because  of  the  degree 
of  freedom  of  cross  section,  are  consequent  on  the  previous  ones;  in  this 
manner  torques  distributed  on  z  will  be  associated  to  a  transversal  elemen¬ 
tary  load,  and  bending  couples  and  bimoaents  will  be  associated  to  axial 
actions.  The  determination  of  inertia  forces  will  complete  the  calculation 
of  the  actions  caused  by  imposed  displacements. 

Such  procedure  is  slstema tical ly  used  for  the  flexural,  torsional 
and  extensional  oscillations  and  permits  to  formulate  the  general  equations 
expressing  the  motion  of  the  thin  walled  berm  in  general  as  well  as  taking 
into  account  the  unstablllzlng  effects  of  stresses  (corresponding  to  the 
basic  position  of  the  beam)  and  of  tne  surface  loads. 


2,3)  Effects  due  to  flexural  motion 


Let  us  consider  the  flexural  deformation.  It  is  characterized 


(fig.  11)  by  displacement  components: 

u  (z,  t)  ;  v  (z,  t)  (2*.  1) 

of  the  line  of  shear  centers  0;  and,  for  the  rotation  of  sections  around 
axes  x  and  y,  by  the  displacement  component  along  axis  z. 


t)  =-  ( 


=  -(«*♦ 


c?v  . 


(2,2) 


We  consider,  above  all,  the  unstabilizing  effects  due  to  stresses 
and  we  calculate,  along  axis  x  of  fixed  system  Cxyz,  the  components 
dfx  of  the  elementary  forces  ac.ing  on  the  elementary  buckled  stripe  dA  dz 
of  the  beam  pertaining  to  two  cross  sections  at  the  distance  dz.  With 
reference  to  fig.  12,  representing  ihe  projection  of  dA  dz  on  the  plane 
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Fig.  11 


x  z,  we  have: 

JlM  < ' * SJ ♦  (*♦  s»' * |f 4l) 

since  the  load  Q°z  (z,  s)  ds  dz  applied  on  the  element  does  not  give  any 
conponent  along  x  and  keeps  the  direction  of  axis  x. 


Developing  equations  (2,3)  we  obtain,  neglecting  quantities  of 
higher  order  than  the  first: 


rig.  is 
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Developing  equation  (2,6)  we  obtain,  for  unit  of  length: 

In  the  sane  manner  as  i< c  equation  (2,6),  we  have  the  moment: 


(2,9) 


</A}x=-S 


tending  to  turn  over  the  element  dA  dz  around  axis  x;  for  unit  of  length, 
we  have  the  elementary  distribution  moment: 

c/'fcr-  -  &T 

*  **  3Z 

and,  for  the  whole  section: 

d\r 


(2,10) 


(2, 10’) 


-t-  * 
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(2,11) 


Integrating  equations  (2,5)  and  (2,9)  cn  the  transversal  area, 
we  have  the  new  distributed  actions  due  to  the  fact  that,  in  buckled 
condition,  normal  stresses  lean  forward  forming  variable  angles  in 
regards  to  the  original  direction  of  z  axis. 


px 


Thus  we  have : 

Equations  (2,5)  and  (2,9)  give  the  transversal  load  due  to  the 


=  1 


13  (*■  dS. 


Jti) 


(2,12) 


flexural  buckling  of  the  elementary  stripe  dA  dz ;  consequently,  we  have 
the  following  twisting  elementary  moment  distributed  along  z: 

***  *  k  fc  ft  [sifi 1 


(2,13) 


using  the  symbols  of  fig.  13  which  shows  as  positive  the  twisting  moment 
(or  the  angle  £  )  if  its  direction  of  rotation  is  the  same  bringing  axii. 
x  on  axis  y. 

Integrating  on  the  whole  cross  section  A  we  obtain: 


(2,14) 


Equations  (2,8)  and  (2,10)  are  always  valid  if  the  loads  Q°z  (z,s) 
keep  the  same  direction  of  axis  z  of  fixed  system  Cxyz.  Let  us  consider 
now  the  effects  of  shearing  stresses  and  actinK  on  the  transversal 

sides  dA  of  the  elementary  buckled  stripe. 

With  reference  to  fig.  14,  showing  the  projection  of  element  dA  dz 
on  plane  x,  z,  we  calculate  the  components  along  z  of  elementary  forces 
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Fig.  14 


acting  on  the  buckled  stripe. 

Therefore  we  have: 


'4  = 


(2,15) 


where  the  effect  of  surface  loads  Q°x  (z,  s)  is  not  present,  since  such 
loads  remain  parallel  to  the  axis  x  of  fixed  system  Cxyz. 

In  conclusion,  for  unit  of  length,  along  z  we  have  the  following 


elementary  axial  load: 

In  the  same  manner,  considering  the  projection  of  buckled  stripe 
dA  dz  on  the  plane  zy  to  calculate  the  effect  of  oblique  in  regards 

of  fixed  axis  y,  we  have  (fig.  15) : 


(2,16) 


**  K  *y  #2  J  ^ 

where  is  not  present  the  effect  of  conservative  loads  Q°y  (z,  s)  which 


(2,17) 


keep  their  direction  along  y. 

Developing  equation  (2,17) .  we 

JO  __  <?  /•*.  01/" 


^  cy  dz.  ' 


(2.18) 
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representing  the  new  action  along  z  due  to  the  different  slopes  of  stresses 
in  the  buckled  state  of  element  <±A  dz;  for  unit  of  length,  from 
equation  (2,18)  we  have: 


Jp-.  ~(t 


(2,19) 


Equations  (2,16)  and  (2,19)  refer  to  the  elementary  area  dA; 
for  the  whole  area  A  of  the  beam  cross  section,  we  nave  the  following 
new  axial  action: 


4  ■  -1M 


<?<C  * 

**  4  *y  sz 


')4V 


(2,20) 


To  equations  (2,16)  and  (2,19)  are  associated  some  distributed 
bending  moments,  since  they  act  at  distance  x  and  y  from  the  axis  of  the 


centroid;  therefore,  we  have  for  ( 
moments  drn^  and  d/3y 


«< 


the  elementary  distributed 


Y  dZ*  **02  ' 


(2,21) 


and  for 


V 

*Y 


9zK  J 

Integrating  on  area  A  we  finally  have: 


rf#r-.f.(r  &*</«) 


(2 , 22) 


(2,23) 


(2,24) 


(2,25) 


furthermore,  equation  (2,16)  gives  the  blooment  variation: 

t?z  J4  3  2^  **■  Jz  *y  **  ' 

Equations  (2,7),  (2,11),  (2.12),  (2,14),  (2,20),  (2,23), 

(2,24),  and  (2,25)  represent  the  new  actions  due  to  the  variable  slopes 
of  normal  and  shearing  stresses  in  the  buckled  stat",  but  it  is  essential 
to  notice  that  in  such  condition  the  forces  acting  on  the  element  have 
a  different  position  if  compared  to  the  fixed  axis  Cxyz.  Obviously,  this 
changes  the  stresses  field  in  the  body;  in  order  to  calculate  this  effect 
it  will  be  sufficient  to  refer  to  the  elementary  stripe  jA  dz  and  consider 
the  moments,  relative  to  the  forces  acting  on  two  sides  dA  as  well  as 
those  actir.-j  on  lateral  surface  of  dA  dz,  due  to  the  displacement  of 
such  forces  frcm  basic  position  to  the  displaced  one. 

We  begin  by  considering  the  effects  of  the  displacement  of 
elementary  shearing  forces  and  ,  iiistributed  on  A, 

and  of  surface  forces  Q®  (z,  s)  and  Q°  (z,  s) . 

A  y 


t: 

V 

tr  cffi 

£X 


s  ^  Q° 

6-<jV5  r 

Z 


-C 

- 1 

tf,  i 


ir! 


7 

1  dZ 


.  Sx'jj, 


y  w 


Fig.  16 


With  reference  to  fig.  16  we  have  for  the  elementary  stripe 
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dA  dx  tne  following  change  of  twisting  nomem: 


-  Ji**  ,J"{  Qfei)  diJiK  -  -j~  + 

*  j|  (t,y (?*(*/ *)<i\dzv-  *  Qy(x,6)d$d2  <c 

Integrating  on  area  A  and  on  center  line  s  of  the  cross  section, 
we  obtain  for  unit  of  length: 


^  V j if, *)  «  o'i  (2>. 

Also  because  of  the  rotation  of  cross  sections  around  axes 

x  and  y,  the  elementary  internal  forces  <•*  ■  and  the 

external  surface  loads  Q°  (2,  3)  aad  Q°  (2,  3)  move  their  points  of 

*  y 


application  of  quantity: 


a*.  ax  ' 


(2,28) 


With  the  same  procedure  previously  used,  we  obtain  the  following 
distributed  elementary  bending  couples: 

**  mk  Cv  c  + 


(2.30) 

(2.31) 

(2.32) 


* 


The  calculation  of  effects  due  to  the  displacement  of  forces 
directed  along  axes  x  and  y  is  now  complete;  now  we  want  to  consider  the 
effects  of  displacements  of  normal  forces  rM  ,  distributed  on  A, 
a”d  of  surface  forces  Q°z  (z,  s) . 

With  reference  to  fig.  16,  we  have  the  following  elementary 
displacement  couples: 


K  2  S2.  * 


(2,33) 


AH  =  S-Art^dz  *  ~‘u  dfidx*(f  (*,*)«  di 

y  *  <?*  ee  * 


(2,34) 


integrating  on  area  A  and  on  line  s,  we  obtain  for  unit  of  length: 

J*  J*  K 

^  _  ( (r  ^  ■<  f  ^ 

y  K 

Now  we  calculate  the  corresponding  inertia  reactions. 

Being  the  mass  for  unit  of  volume  of  the  thin  beam, 

for  the  elementary  mass  dA  dz  the  following  forces  correspond  to  the 
displacements  (2, 1)  and  (2,2): 

<tL=  -udddz  & 

*  1  dtl 

..  J l J  J c 


obviously,  to  (2,37)  correspond  the  distributed  couples: 

=-IjM  (mF*X  ^  *  /•  d^y=rJA{M^ 


(2,35) 


(2,36) 


(2,37) 


(2,38) 
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Integrating  on  the  whole  section  we  obtain  the  following  actions 


for  unit  of  length: 


Py r* 


-  o 


(2,39) 


40) 


On  the  contrary,  the  biaoment  which  seems  to  develop  from  (2,37) 
equals  zero;  in  fact,  we  have:  ^ 

1 .  u  If  G^k-  jc  *  “  y)  -  C? 

■gj’  *  t*  '  (2,41) 


because  the  sectorial  coordinate  t>j  is  orthogonal  to  the  coordinates 
x  and  y. 

2,4)  Effects  due  to  torsional  oscillation 

Let  us  consider  the  torsional  buckling  shown  in  fig.  17. 

Since  the  cross  sections  rotate  around  the  shear  center  axis, 
every  element  of  the  area  dA  moves  along  x  and  y  as  follows: 


u  (z,  t)  =  (yQ  -  y)  £ 


(z,  t)  =  -  (xQ  -  x)  P 


and  moves  along  z,  because  of  the  warping 

t?  &> 

w  (z,  t)  =  -  W 

as  it  results  from  the  sectorial  areas  theory. 


(2,42) 


(2,*3) 
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The  angle  ^  ,  together  with  the  twisting  moment  Mz,  is 
therefore  considered  positive  if  it  brings  x  on  y,  being  z  downward. 

Considering  above  all  the  effect  of  stresses  gT  we  calcu- 
.ate  the  components  along  x  of  elementary  forces  acting  on  the  buckled 
form  of  the  elementary  stripe  dA  dz  (fig.  18). 


As  in  the  case  of  flexural  motion,  we  obtain  the  transversal 
action  relative  to  the  buckled  element  dA  dz: 

df=  -6;  (6^.  (  $  *  f £**) 


(2,44) 


since,  also  in  this  case,  the  load  Q°^  (z,  s)  keeps  its  direction.  From 
equation  (2,44)  we  obtain,  with  reference  to  the  unit  of  length,  the  fol¬ 
lowing  elementary  transversal  forces: 


(2,45) 


In  the  same  manner,  considering  the  projection  of  dA  dz  on 
plane  y  z  we  obtain  (fig.  19): 


(2,46) 
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Integrating  on  A,  we  obtain: 

'ythfcv*)%l',A 


(2,47) 


(2,48) 


which  represent  the  transversal  loads  developing  on  the  thin  beam  slightly 
twisted  due  to  the  different  slope  of  stresses  acting  along  the 

fibers.  The  loads  Q°z  (z,  s)  also  in  this  case  remain  parallel  to  axis 
z  of  fixed  system  Cxyz. 

Furthermore,  as  in  section  2,2),  for  equations  (2,45)  and  (2,46) 
we  have  elementary  turnover  moments  d  and  d  My  due  to  the  components 
along  y  and  x  of  elementary  forces  *  aW  ;  they  are: 

a  (2,49) 


(2,50) 


Integrating  equations  (2,49)  end  (2  50)  on  A  we  obtain,  for  unit  of  length: 

i  *<(*.-*•>  If  O.M 


*,  (y.-V 

*  result  of  equations  (2,45)  and  (2,46)  we  obtain  the 
elom^ntary  twisting  moment  distributed  as  follows: 


(2.51) 

(2.52) 


(2,53) 


duo  to  the  fact  that  dp^  and  dpy  act  at  distances  (yQ  -  y)  and  (xQ  -  x) 
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» 


* 


fro*  ttM  axis  of  shear  cantar;  integrating  me  obtain: 


(2,54), 


Lea  u a  conaider  now  the  unstabilizing  affects  due  to  ahaarinf 
strasaaa  and  which,  because  of  toraional  buckling,  produce 

componenta  along  the  directions  of  azea  x,  y  and  x  of  fixed  system  Cxyz. 
Since  the  fibers  of  the  thin  bean  bend  because  of  (D{x ),  as  we  did  for 
tho  flexural  notion,  we  calculate  the  component*  along  x  of  elementary 
forces  and  ,  distributed  on  A  and  variable  along  the 

buckled  fiber. 


forces 


did  for 


Kroa  the  projection  of  element  dA  dz  on  planes  zx  and  zy,  we 
obtain,  in  accordance  with  figs.  14  and  15: 

Integrating  on  A,  for  unit  of  leugth,  we  obtain  the  axial 


distributed  load 


ui«u  iuia  j 

4  "■]»/«  Y.'Y)-^*.-*>]d*\ 

Equation  (2,56)  furnishes  the  distributed  bending  moment?: 

j  g&h-v-wsH* 


(2,57) 


(2,58) 


and  the  hi mo went  variation: 


§}--j£lU  o-  %  '■  d”\ 


(2,59) 


Las  us  calculate  now  the  unatabilizing  notions  corresponding 
to  components  of  *^4 A  and  ,  distributed  on  A,  dn  axes  x  and  y. 

We  cbnsider  tho  elementary  stripe  dA  dz  of  fig.  20;  It  will 
be  stressed  on  two  sides  dA,  of  abscissa  z  +  dr  and  z  respectively,  by 


Fig.  20 


elsnentarj  forces 


°<)*  Vt* *>**  '  ( V  J*H' :sW-  % ^  /  $r  ^ 

ben F  in  regards  of  fixed  axes  x  and  y  of  »  ■*lere ^ J  and 

,  boing  j-  and  £  tbe  direction,  of  aovefcla  uu  x  and  y. 

The  loads  Q“x  (z,  b)  dedr  end  Q°y  (t ,  a)  dad*  thl  act  on  tbe  aide 
surface  Aa  dz  of  tbe  stripe  and  will  «*-»?  the  directions  of  fixed 
exaa  x  and  j . 

By  c  niculating  th#  components  on  x  and  y  of  all  forces  seeing 
on  dA  «£z,  wo  obtain  the  transversal  loads.  Thus  we  have: 

*<£•(**/  *£*^4*  (}<%-<*/ . 

^ C^v fcr^  4$  fut  jp  x 


"  $]?? A lr  <M  2£dt*<f(%,)Jidt 

Si  fix  f  /  *f  •* 


(2,80) 


For  the  equilibriua  in  the  basic  position,  we  have; 

f  <?°  d*  -T  O 

■* 

and  equation  (2f60)  ia  simplified  as  follows: 


<4  = 


(2,61) 


representing  the  new  transversal  forces,  distributed  on  dz  and  directed 
along  x,  due  to  the  new  slope  of  shearing  stresses  after  buckling;  for 
unit  of  length,  we  have: 


being  c  . 

In  the  mm  Banner;  performing  a  projection  along  y,  we  obtain 
the  elementary  transversal  forces  directed  aloLg  y  and  relative  to  the 
length  dz  of  element  dA  dz; 


being,  for  the  equilibrium  along  the  direction  j  in  the  basic  condition 


*  Qy  =  © 

Sr  Jf  ' 


For  unit  of  length,  we  have: 


(  ^  f 


(2,63) 


Integrating  aquations  (2,62)  and 
section,  we  obtain  the  transversal  leads: 

l*-Jk^***> 

4 


(2,6r)  on  the  area  A  of  cross 


12,66) 


rY  » 

Furthermore,  the  components  aljng  x  and  y  of  and 

*1 

produce  the  following  turnove:.  couples  around  ares  x  and  y 


(2,67) 


(2,68) 


4  hK  »  -  tr  jdQd*. 
4^-  $4*4i 


Integrating  os  A,  n  obtain,  for  unit  of  length: 


<*r  -  L  ^ty  f  d$ 


(2,66) 


(2,70) 


(2,71) 


From  aquation*  (2,62)  and  (2,65)  w*  obtain  th*  distributed 


tJtmitary  touting  moments ; 


(2,73) 


*  dP,  ( y.-y)  -  <ipr  (*.-*)  - 

(2  72) 

^df>)  (y.-y;  ~~ 

from  which,  bp  lntagratlon  on  araa  A,  we  have: 

In  thia  wanner  we  have  calculated  the  new  action*  developing 

•x. 

along  abscissa  z  of  the  thin  walled  beam  In  a  slightly  buckled  fora  due 
to  the  variable  slopes  of  th*  stresses.  For  the  calculation  of  such 
actions  w*  did  not  consider  the  warping  of  the  cross  section,  becauae 
It  causes  only  a  variation  of  generatrices  length  and  not  their  bending; 

4n  the  contrary,  shown  later,  It  will  affect  the  calculation  of 
displacement  moments. 

For  the  calculation  of  such  displacement  moments,  we  consider 
that;  because  of  displacement* 

u  *  (y0  -  y)  ^  ;  v  «  -  -  X)  f 

the  mb  tiing  stresses  end  ^  and  the  surface  loads  Q°x  and 

Q  moved  along  y  and  x  reaoectlvely ;  thus  we  have,  as  in  equation 
(2,27)m  the  elementary  distributed  twisting  moments: 

= £  h,  (**'*>  t  <**]  *£  C%  (Yc~vt  **] 

(2  7i 

*  0%  v  i)  (4- y)  j  <** 


(2,74) 


from  which,  Integrating  on  A,  we  obtain: 


(2,75) 
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*Jt  *  Qr^iyty.-v)]  4  ■ 

L«t  u*  consider  now  the  warping  of  cross  section  of  positive 
sectorial  coordinate;  because  of  this  warping  .which  varies  along  s 
with  ^  ,  the  shearing  stresses  and  aad  th*  *urface  loads 

Q°z  (z,  s)  tsd  Q0^  (z,  s)  wove  in  parallel  with  themselves;  in  tee  save 
wanner  as  for  equations  (2,31)  and  (2,32),  we  have  the  distributed 


bending  couples : 


"  5 

Js> 


(2,78) 


(2,77) 


Nov  we  calculate  -he  effect  of  *41  ajuacement  of  normal  stresses 
6^  and  of  surface  loads  Q°.  (a,  a). 

As  for  equations  (2,33)  and  (2,34)  we  have  the  displacement 
elementary  couples: 


g  o.~*;  W*** 


<///,*  ti&CY.-yjfi*** 

< 

f roe  which,  fcy  Integrating  on  A,  we  obtain: 


(2,78) 


(2,78) 


/ <$  (y.-y>g*.[ f <?> »X/. -yj 4  d* 
i  Z  4  4  (2.1 


Now  we  calci’laL.j  the  Inertia  forces  appearing  during  the 


torsional  motion. 


With  reference  to  oiomentary  mass  dA  da,  we  have  the  elementary 


force 
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<  ty  -  j4  <MJz  (x.'XJ  |p 

<* 

correspondingly,  ve  have  the  elementary  couples: 

^  =  -A  **  •«  f&i 
</*,.-  -/<  4fi<kjyc»*)*it.-y/]  df- 


(2,82) 


(2,83) 


Equations  (2,82)  and  (2,83)  give  the  cross  and  axial  distributed 


forces : 


w* vwu  •  a 

f>  T-HQy 

*k  r  '•  ah 


p*=° 


(2,84) 


and  the  couples : 


=o 


Wx  =  o 


(2,85) 


while  the  components  df2  give  the  bimoment  variation: 

j6  r 

7iz  f4  “  573?*  <2'86> 

being  3^  and  XuS  the  quantities: 

L=[1*.-x)'hy.-y)‘]jq = 4  *ty*0(vY.1);  <2  S7) 

T9  ^ 

polar  moment  of  the  cross  section  in  regards  of  the  shear  center  0  and 
sectorial  moment,  respectively. 


2,5)  Effects  due  to  extensional  oscillation 

Let  us  consider  the  extensional  oscillation. 

Since  the  sections  have  only  displace  tents  w(z)  along 
(fig.  21) ,  in  the  extensional  buckling  we  do  not  hare  variable  slopes 
of  stresses,  in  regards  of  axes  x,  y  and  z,  and  corresponding  unstabi¬ 
lizing  effects  do  not  appear.  On  the  contrary,  we  notice  some  unstabi¬ 
lizing  effects  because,  due  to  the  extensional  displacement  w,  variable 


with  z,  th«  shearing  stressea  "t”  and  and  the  loads  Q°,  (s,  s) 

<Jt  m7 

and  Q °y  (z,  s)  move  In  parallel  with  themselves,  So  w*  have  the  bending 
distributed  couples: 


Tig.  21 


The  inertia  forces  are  given  by:  ^ 

Pj.  *  0  Py  -  0  PE  -  fl*  (2,90), 

■^*^■■^■8-0  (2,91) 
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2, *}  Several 


equations  of  dynamics  o f  thin  walled  beams 


In  paragraph.  2,2),  2,3)  and  2,4)  above  we  calculated  the  actions 
develop!.,  on  the  thin  mailed  beam  during  it-  displacement  caused  by 
flMUral,  torsional  and  extension.!  buckling.  They  are  alv.y.  balanced 
*ith  th*  elastic  reaction,  and  the  loads  directly  applied;  if  we  approxl- 
aw9e  the  curvatures  in  the  planes  *r  and  yz  to  the  curvatures  in  }$ 

**  ji  ’  we  wl11  obtfiln  the  following  differential  equations  system: 

']f(§  -ji  J£  <%**) 

-r*rS  -J  tSU 

■0S X*«*J 

*  -9*  a,  *  JS  (2,92) 

,  -  fa*  -j pj y4i 

1  -  trt/  (X.-XJ] 

>fr;cvxjf  QfayjJM 


!&* 
4 


j<&  5 '  %  T>«]  -[&{§£ kx  (y.-yj -y v-#W 


ef*  *  P-h 


*  <& 
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r»lr.  &*&}  * 


'  J 


+ M  t&'O  -  Ql (r^  4i  +I  Tt*  * Ji 


Above  #y»tes,  together  with  the  boundary  and  initial  condition#, 

furnishes  the  action  of  the  bees  under  (F°)  and  (<&*) ;  these  letters  are 

represented  directly  by  distributed  loads  Q°  (z,  a),  Q°  (z,  s)  and 

*  / 

Q°r  (a,  s)  and  indirectly  by  stresses  component#  -  **  ■ 

corresponding  to  the  basic  equlllbriua  condition. 

As  far  as  stability  is  concerned,  it  la  interesting  to  d  teralne 
the  value  of  the  multiplier  ^  of  (F°)  and  (  £*)  for  which  the  motion 
is  no  acre  limited  jtosuch  case  of  conservative  forcea  and  dislocations 
the  change  from  stability  to  unstability  will  be  expressed  by  the  value 
zero  of  the  motion  frequency. 

Above  system  (2.93)  also  includes  all  the  problems  of  stability 
and  dynamics  of  thin  walled  beam  of  close  section  end  of  the  solid  section 
bean  and  can  be  easily  applied  to  the  various  particular  cases,  expressing 
from  time  to  time  the  applied  loads  and  the  stresses  components 
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